The /c-local Pauli Commuting Hamiltonians Problem is in P 
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Given a Hamiltonian that is a sum of commuting few-body terms, the commuting Hamiltonian 
problem is to determine if there exists a quantum state that is the simultaneous eigenstate of all of 
these terms that minimizes each term individually. This problem is known to be in the complexity 
class quantum Merlin-Arthur, but is widely thought to not be complete for this class. Here we 
show that a limited form of this problem when the individual terms are all made up of tensor 
products of Pauli matrices is efficiently solvable on a classical computer and thus in the complexity 
class P. The problem can be thought of as the classical XOR-SAT problem over a symplectic vector 
space. This class of problems includes instance Hamiltonians whose ground states possess topological 
entanglement, thus showing that such entanglement is not always a barrier for the more general 
problem. 



I. INTRODUCTION 

Finding the energy of the ground and excited states of 
a many-body quantum system is a canonical problem in 
physics. While there are many simple cases where one 
can exactly analytically solve for the energies, in general 
one often needs to resort to computer simulation to get 
even approximate answers to these problems. Even then 
it is often the case that this problem is computationally 
intractable as a function of the size of many-body prob- 
lem. Recently, motivated in large part by the birth of 
quantum computing (whose original motivation was in 
part due to the intractable nature of simulating quan- 
tum systems |l| ) , a new understanding of the computa- 
tional complexity of these quantum many body problems 
has begun to emerge. A seminal result in this literature 
is that of Kitaev [2], Q who showed, roughly, that the 
problem of determining whether a system has a ground 
state energy below or above an energy level (specified to 
polynomial accuracy) is complete for a quantum equiv- 
alent of the computational complexity class NP (non- 
deterministic polynomial time), the class QMA (quantum 
Merlin- Arthur). This means that the problem is both as 
hard as every problem in QMA, and that if it could be 
solved efficiently it would allow for efficient algorithms 
for every problem in QMA. Or, more loosely, that it is 
unlikely that even a quantum computer could help with 
this problem. 

Since Kitaev's pioneering work tremendous progress 
has been achieved in trying to understand the conditions 
under which this problem remains QMA-complete [4|-[llj . 
For example if the Hamiltonian of the many-body system 
is made up of a sum of commuting two-qubit or three- 
qubit interactions, then the problem is NP-complctc in- 
stead of QMA-complete [a, [llj. The class of problems 
where the input Hamiltonians is made up of a sum of 
commuting terms has drawn special attention because 
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these Hamiltonians can exhibit non-trivial entanglement 
(topological entanglement [12| , for example) and yet ac- 
cording to physics lore should be thought of as essentially 
classical [a, [UJ, \}A\. Here we add another small result 
to this pantheon of results by showing that a quantum 
many-body system whose Hamiltonian is a sum of com- 
muting terms that are made up of tensor products of 
Pauli matrices is solvable in polynomial time, i.e. is in 
the complexity class P. While this is perhaps not sur- 
prising to those familiar with the theory of stabilizer er- 
ror correcting codes [13l |. the Hamiltonians we consider 
include systems that exhibit topological order, such as 
Kitaev's toric code [Tj], thus perhaps shedding light on 
the role that topological order plays in the difficulty of 
the commuting Hamiltonian problem. Our result can be 
thought of as instances of XOR-SAT over a symplectic 
vector space. 



II. THE LOCAL HAMILTONIAN PROBLEM 

Here we briefly review the local Hamiltonian problem. 
We begin by defining the problem: 

Problem 1 The fc-local (qudit) Hamiltonian problem 
Given: A set of Hermitian positive semi-definite oper- 
ators, {Hi, 1 < i < r} acting on a Hilbert space of n au- 
dits (fixed dimension d) each of which is bounded norm, 
\\Hi\\ < c ; f or some constant c. Each term Hi acts non- 
trivially on at most k of the n qudits and is specified to 
accuracy polynomial in n. Additionally we are given two 
numbers, b > a that are separated by an inverse polyno- 
mial gap, b — a > l/poly(n). We are promised that the 
smallest eigenvalue of H = ^■ 1 Hi is either smaller 
than a or greater than b. 

Decide: Determine whether the smallest eigenvalue of 
H = ^2 i= i Hi is smaller than a or greater than b. 

The fc-local Hamiltonian problem is meant, in part, to 
capture a natural class of problems that are encountered 
by physicists. To this end the problem can be naturally 
thought of as attempting to find the energy of the ground 



state of a quantum-many body system. Note however 
that the problem does not require actually producing an 
efficient description of the ground state. 

The fc-local Hamiltonian was introduced by Kitaev 0, 
|3j who showed that, with certain restrictions on the set of 
Hamiltonians, this problem is complete for the complex- 
ity class QMA. Problems in QMA are, roughly problems 
for which there exists a (quantum) proof such that if an 
instance is in the language then it can be efficiently ver- 
ified on a quantum computer. Thus Kitaev showed that 
for every "yes" instance of the fc-local Hamiltonian prob- 
lem there is an quantum state on poly(n) qudits which 
can be used to verify that this instance is a "yes" in- 
stance. Conversely Kitaev also showed that every prob- 
lem that has an efficient quantum proof that can be ver- 
ified in polynomial time on a quantum computer can be 
reduced into an instance of the fc-local Hamiltonian prob- 
lem. 

Kitaev's original result on the QMA-completeness of 
the fc-local Hamiltonian problem placed certain restric- 
tions on the Hamiltonian. In particular his proof only 
held for qubit Hamiltonians with k > 5. Subsequently 
this work was improved to fc = 3 (qubits) [7|] and fc = 2 
(qubits) 0. In another direction it was shown that the 
fc = 2 result holds even for Hamiltonians whose interac- 
tions graphs are planar [9( , and even for quantum system 
in one spatial dimension when using qudits instead of 
qubits {d = 12 in jj|). Thus a wide swath of Hamil- 
tonians result in fc-local Hamiltonian problems that are 
QMA-complcte. 

In a different direction one can also consider restric- 
tions on the class of Hamiltonians being considered which 
reduce the complexity of the problem to (presumably) 
weaker complexity classes like NP and P. For example if 
the Hamiltonian is diagonal in the computational basis of 
each qubit, then a seminal result of Barahona [16[ shows 
that the problem is not QMA-complete, but instead is 
NP-complete. Given this fact it is interesting to con- 
sider Hamiltonians that arc made of sums of commuting 
fc-local terms. These can be simultaneously diagonal- 
ized so that in some basis they effectively behave like the 
classical problem, and yet because the diagonalization is 
non-trivial it is not clear that all of these problems are 
NP-complete. This was investigated by Bravyi and Vya- 
lyi Q who showed that if the Hamiltonian is 2-local then 
the problem is NP-complete. Recently Aharonov and El- 
dar [Tl| showed that this holds for 3-local commuting 
Hamiltonians as well, for both qubits and qutrits. This 
is especially nice because it shows that these systems do 
not possess topological order and in effect possess only 
local entanglement. Motivated by this result in this pa- 
per we investigate a class of Hamiltonians that does pos- 
sess topological order, those that are like Kitaev's toric 
code [14[ and made up of a sum of commuting Pauli op- 
erators. In contrast to all of the problems above we here 
show that this problem is actually efficiently solvable in 
polynomial time. 



III. THE COMMUTING HAMILTONIAN 
PROBLEM 



Having reviewed the fc-local Hamiltonian problem, let 
us define the promise problem, fc-local (qubit) Commut- 
ing Hamiltonian Problem: 

Problem 2 fc-local (qubit) Commuting Hamiltonian 
Given: A set of commuting k-local projectors on n qubits, 
{Hj},l < j < r = poly(n) each of whose entries are 
specified using poly(n) bits of precision, and each of 
which pairwise commutes [Hj.Hk] = for all j,k. We 
are promised that either 

(a) There exists a zero energy eigenstate of H = 
Y^ij=iHj> *- e - ^W) suc h that H\ip) = 0. 

(b) There is no zero energy eigenstate and the lowest en- 
ergy of eigenstates is at least 1. i.e. \/\ip) there exists 
a j such that (ip\Hj\ip) > 1 (this could be inverse 
polynomial, but in the case of commuting Hamiltoni- 
ans the lowest energy, if it is not 0, can be made into 
an equivalent problem with the lowest energy > 1.) 

Decide: Whether the given instance obeys condition (a) 
above. 

The commuting Hamiltonian problem is known to be 
NP-hard but is not known to be NP-complctc. A com- 
monly used counter-example to the idea that commut- 
ing Hamiltonian problem is NP-complete are Hamilto- 
nians such as Kitaev's toric code Hamiltonian. These 
Hamiltonians are made up of commuting fc-local terms 
yet possess a large amount of entanglement. In partic- 
ular these Hamiltonians possess a global order, topolog- 
ical order, and hence cannot be prepared by a constant 
depth quantum circuit. However, here we will show that 
this intuition is wrong and that the commuting Hamilto- 
nian problem where all Hamiltonians that are made up of 
commuting projectors onto eigenspaces of fc-qubit Pauli 
operators belongs, like the models considered by Kitaev, 
in fact, to the class P. 

Define the following problem, of which Hamiltonians 
like Kitaev's toric code are a particular example: 

Problem 3 fc-local Pauli Commuting Hamiltonian 
Given: A set of commuting k-local projectors on n qubits, 
{Hj}, 1 < j < r = poly(n) each of which pairwise com- 
mutes [Hj, Hk] = for all j, k and for which each Hj is 
a projector onto the 1 eigenspace of a k-local Pauli oper- 
ator, Hi = o(I ~ Si) where Si is a k-local Pauli operator 
with Sf = I . We are promised that either 

(a) There exists a zero energy eigenstate of H = 
y] i — 1 Hj, i.e. 3\tp) such that H\ip) = 0. 

(b) There is no zero energy eigenstate and the lowest en- 
ergy of eigenstates is at least 1. i.e. W\ip), (ij>\H\i()) > 



Decide: Whether the given instance obeys condition (a) 

above. 

Claim: k-local Pauli Commuting Hamiltonian problem 

is in P. 

Recall that the single qubit Pauli operators are / = 



X = 



Y = 











and Z = 




-1 



and that Pauli operator on n qubits is a tensor product 
of these operators (with a possible phase a multiple of i). 

The proof of this claim will now occupy the rest of 
this section. The main intuition is that this problem is 
related to a question about stabilizer codes, a question 
for which a sufficient Clifford circuit can transform this 
problem into an XOR-SAT problem, which is efficiently 
solvable on a classical computer. Further this Clifford 
circuit can be efficiently computed and simulated. Thus 
the essential tool is the Gottesman-Knill theorem [17| ■ 

As a first step we recast this problem in a slightly more 
palpable manner for those who know stabilizer codes. In 
particular if we write Hi = Ml— Si), where Si is a fc-local 
Pauli operator that squares to identity Sf = I, then the 
two conditions become 

(a) 3\i>) such that Vi, Si\ip) = \4>). 

(b) V\ip) that are common eigenstates of the Si's, there 
exists at least one i such that Si\ip) = — \tp) 

The problem now is given Si that are fc-local Pauli op- 
erators, decide whether there is a common +1 eigenstate 
of these operators. Here we will show how this can be 
done in polynomial time on a classical computer. To be 
precise, define the 

Problem 4 fe-local Pauli Commuting Hamiltonian (Sta- 
bilizer Version) 

Given: A set of commuting k-local Pauli operators, {Sj}, 
1 < j < r = poly(ji) which all square to identity. We are 
promised that (actually it is always true that) 

(a) 3\ip) such that \/i, Si\tjj} = \tp). 

(b) y\ip) that are common eigenstates of the Si 's, there 
exists at least one i such that Si\ij)) = —\ip) 

Decide: Whether the given instance obeys condition (a) 
above 

At this point we should note that we can further drop 
the condition of fc-locality as so we will drop this from 
here on out. 

We begin by reviewing some basic facts. Every Pauli 
operator that squares to identity can be written as 
(-l) p ®" =1 P{xi, Zi) where, p e {0, 1} and 



P{Xi,Zi) 



I if Xi = Zi — 
X if Xi = 1 , Zi = 
Y if Xi = Zi = 1 



(1) 



Z if Xi 



0,Zi 



1 



Thus we can associate to every such operator a In + 1 



our r Si operators we can denote these operators us- 
ing r 2n + 1 binary vectors, which we will denote 
(Pi,Xi,i,---,%i,n,Zi,i,---,Zi,n)- We can assemble these 
into a r by In + 1 matrix, which we call M. 

Now notice the following stability of our problem. Sup- 
pose >S = {Sj\l < j < r} is an instance of our prob- 
lem. Let S' denote the set formed by replacing the 
kth element of S by SjSk where j ^ k. We claim 
that S satisfy (a) iff S' satisfy (a), and similarly for 
(b). To see this note that if S satisfy (a), then 3\ip) 
such that \/i Si\tp) = |i/>), and hence for this same \ip), 
Vi 7^ k,Si\ip) = \i/j), while SjSk\ip) = \ip), since we can 
apply each Sj , Sk in turn. Similarly if S' satisfy (a) , then 
3|V>) such that Si\i>) = \ip), i ^ k, and SjS k \i)) = \ip), 
and hence for this same \tp) Si\ip) = \ip) for all i, with 
the notable case that Sj(SjSk)\i>) = Sk\ip) — \i>)- 

For the no instances we proceed similarly. For conve- 
nience let's denote S' by S' = {Si\l < i < r and i ^ 
k} U {SjSk}, where j ^ k. If S satisfies (b), then V| V 7 ) 
that are common eigenstates of Si's, there exists at least 
one Si & S such that Si\ip) = — \ip). If i ^ k, then we 
are done, since there is at least this same Si in S' such 
that Si\i/>) = -\tp). If i = k, that is S k \ip) = -\ip), then 
there are two possibilities. One possibility is that Sk is 
the only one in S such that Sk\ip) = —\ip) and implicitly 
for all other Sj £ S, Sj\ip) = \ip). Then we would have 
SjS k G S' such that SjS k \ip) = -Sj\ip) = -\tp). The 
other possibility is that there are other Sj's such that 
Sj\ip) — —\ip)- Then we have at least those Sj's in S', 
no matter whether SjSk\ip) = — \ip) or SjSk\ip) = \ip)- 
Hence if S satisfies (b), then S' satisfies (b). In the other 
direction, if S' satisfies (b), there are also two possibil- 
ities. If there is any Si 6 5', where i ^ k, such that 
Si\ip) = —\ip) then we are done. Otherwise if SjSk is 
the only one in S' such that SjSk\ip) = — \ip) and implic- 
itly for all other Sj E S' , Sj\ip) = \tp). Then we have 
Sfc|'0) = — \i J ) m S. Hence we have shown that if S' 
satisfies (b), then S satisfies (b). 

Thus we have shown that an instance S can be con- 
verted into a completely equivalent instance, S 1 by re- 
placing an element of S with a product of that element 
and another element from S. Suppose that M is the 
matrix describing the S as above and M' is the matrix 
describing S' . Ignoring, for the moment, the first col- 
umn of M (corresponding to the phase factors), M' is 
simply the matrix for M with the fcth row replaced by 
the bit-wise addition modulo 2 of row j and row k of 
AI. The phase factor can be calculated as follows. First 
note that only the phase factor for the kth row of M' is 
(potentially) different from that of M . Define 



g(xi,zi,x 2 ,z 2 ) 



if x\ = Z\ = 

z 2 — x 2 if x\ = z\ = 1 

z 2 (2x 2 - 1) if X\ = 1, z\ = 

k £2(1 - 2z 2 ) if X\ = 0, Z\ = 1 



(2) 



vector of binary numbers (p, X\, 



, x n , Z\ , 



z n ). For 



Then if 2r 3 + 2r k + Y^=i 9(xj,i, Zj,i, x k ,i,Zk,i) = mod 4 
set the first bit of the kth row of M' to otherwise set the 
first bit to 1 (the sum will only ever give 2 mod 4.) While 



this phase factor is a bit messy the only real important 
factor is that in going from S to S' we can efficiently 
calculate the new phase factor of M' from M. 

Next note that our problem is invariant under unitary 
conjugation of the {Sj} for a particular instance. That is 
if we let S = {Sj\l <j< r}, and S' = {USjW, 1 < j < 
r}, then S satisfies (a) iff S' satisfies (a) and similarly 
for (b). This follows simply from the fact that unitary 
conjugation corresponds to a change of basis and both 
statements of (a) and (b) are basis-agnostic. Further if 
the U used is a Clifford group unitary, then the new set 
of operators S' will all be elements of the Pauli group and 
thus in this case we obtain an equivalent Pauli Commut- 
ing Hamiltonians problem. Suppose that v is a 2n + 1 
dimensional binary vector corresponding to a Pauli group 
element. Then all elements of the Clifford group can be 
written as a linear transform over Z2 that acts on the 
last 2n elements of v along with a transformation on the 
first bit (the global) phase that is a function of the el- 
ements of v (this transform is not linear, but is easy to 
calculate.) Not all linear transforms can be implemented 
this fashion since the structure of the Pauli group must 
be preserved. 

We will use in particular three transforms, the 
controlled- not C'x, the Hadamard H, and a square root 
of Z gate S = -jx{I — iZ). The Hadamard on qubit i 
simply has the effect of swapping the X and Z portions 
of the vector corresponding to the ith qubit, and changes 
the global phase bit r to r' = r © XiZi where Xi is the X 
component of the ith qubit and z, is the Z component of 
the ith qubit. The Cx gate, acting from the ith qubit to 
the jth qubit has the effect of replacing the jth qubit's 
X component of the vector to the sum of this component 
and the X component of the ith qubit, or Xj — > x% © Xj, 
while replacing the ith qubit's Z component of the vector 
to the sum of this component and the Z component of 
the jth qubit, or z, — >• Zj © Zj. The global phase is up- 
dated as a function of these values in a simple fashion, it 
is the sum of the old global phase plus the product of the 
X component of the ith qubit Xi times the Z component 
of the jth qubit Zj and the X component of the jth qubit 
Xj plus the Z component of the ith qubit Zi plus one, or 
r' = r © XiZj{xj © Zi © 1). Anyway it is easy to calculate 
function of the involved vector components. The Cx, ig- 
noring the global phase, thus has the effect of performing 
two simultaneous row additions, which row being added 
to which depends on which qubit is the target and which 
is the control. The final gate, a square root of Z gate, 
transforms Pauli operators as SXS' = Y, SYS 1 * = —X. 
and SZS^ = Z. Thus if we apply this gate to the ith 
qubit it has the effect of changing the Z component of 
the ith qubit into the addition of the Z component of 
the ith qubit and the X component of the ith qubit, 
or Zi — ^ Xi © Zi. The global phase bit r is updated to 
r' = r @XiZi. 

Now we are ready to describe the algorithm for solving 
the Pauli Commuting Hamiltonian problem, specified by 
a given M matrix. The tools that we have available to 



transform this M to an equivalent problem are the ability 
to add rows, to perform the simultaneous column addi- 
tions corresponding to the C'x operator, the ability to 
swap the X and Z components, to permute both the X 
and Z components simultaneously using a permutation 
gate, and finally the ability to permute the rows of the 
matrix. All of these can be done while efficiently updat- 
ing the global phase bits for the matrix. We will proceed 
by showing that we can use these operations to put M 
into a form from which, either we have discovered along 
the way that the instance cannot satisfy (a) or there ex- 
ists a \ip) satisfying (a). 

Let us denote the M matrix by 



M= [R\\A\\B] 



(3) 



where R is the r by 1 matrix corresponding to the global 
phase bits, and A and B are the r by n matrices corre- 
sponding to the X and Z components of these operators 
respectively. We now proceed as follows: by using per- 
mutations, Cx and row sum, each of which preserves 
whether an instance satisfies (a) or not, we can perform 
Gaussian elimination on the matrix A. When we do this, 
we will end up with a matrix of the form 



Mi 



' Ri 


I 


Bi d" 


i?2 





B2 C2 



(4) 



where / is a k by k identity matrix, with k equal to 
the rank of A. If, during this elimination, we ever ob- 
tain a row that has a global phase bit of 1 and all zeros 
in the remaining rows, then we can immediately answer 
that the instance satisfies (b), because this corresponds 
to having the element —I in the set of commuting Pauli 
operators for this instance, and it is trivial that (a) can 
never be satisfied in this case. Using the fact that the 
operators corresponding to the rows must all commute 
we can further see that it must be the case that B2 = 0: 



Mi = 



' Ri 


I 


Bi C x ~ 


R 2 





c 2 



(5) 



Furthermore, for the same reason, B\ must be a symmet- 
ric matrix which could be further transformed into the 
form NN T , where TV is an invertible matrix UM, 



Mi 



Ri 

R 2 



I 




NN T d 

c 2 



(6) 



If we apply Cx gate according to the form of N, we would 
obtain the new matrix 



Mi 



~Ri 


N 


A^ C x ' 


R 2 





C 2 



(7) 



Now if we apply S gates on the first k qubits, then the 
second N matrix on the first block row would become 
zero. And if we perform Gaussian elimination again on 
the first N matrix, we would transform it into I and get 
matrix M2, 



M,= 



Rn 



I 




d 

c 2 



(8) 



If wc now apply the Hadamard gate to the last n 
qubits this transforms M 2 into 



Mo = 



R > 



I Cx 
C 2 








(9) 



We can now repeat the Gaussian elimination on the X 
component of this matrix, terminating if we ever obtain 
a row or more corresponding to — /. This will produce 



Ma 



R 3 

o 



/ 









(10) 



where / is now a fc' by fc' identity matrix, and we have 
used the fact that the global phase of the lower compo- 
nent must be zero or else we would have terminated upon 
attaining a —I row. Thus we have shown that it is pos- 
sible to transform an instance of the problem given by 
M into an new equivalent instance given by M& where 
we terminate, answering in the negative for the problem, 
if we ever obtained a —I row. Finally we see that this 
problem is nothing more than the commuting Hamilto- 
nian problem for a bunch of single qubit ±X operators. 
This problem always has a solution to condition (a) given 
by an appropriate tensor product of computational basis 
states such that each ±X acts as trivially on this ba- 
sis state. In other words we will always answer yes if we 
make it all the way to the end of the reduction to M4 and 
only answer no if somewhere along the line wc obtained 
a —I row. 

Each of the operations described above can clearly be 
done in polynomial time [la , Il9j . This in turn completes 
the proof that the fc-local Pauli Commuting Hamiltonian 
problem is in P. 



nian is in NP (and hence NP-complcte). Here we re- 
stricted ourselves to the case that when the Hamilto- 
nian is a sum of commuting terms which are made up 
of purely tensor product of Pauli matrices, and proved 
that the fc-local Pauli Commuting Hamiltonian problem 
is actually in P. The fundamental reason for this result 
is that commuting Pauli operators have eigenstates that 
are describablc within the stabilizer formalism which it- 
self is amenable to polynomial time classical manipula- 
tions. This means in particular that for any fc-local Pauli 
Commuting Hamiltonian problem, we could transform 
the individual terms into some basis where the operators 
are all simultaneously diagonalized, and in this basis the 
problem is then a simple XOR-SAT instance. In compari- 
son to other classes of commuting Hamiltonians that have 
been considered this problem ended up in P because not 
only did the eigenstates have an efficient description in 
terms of stabilizers, but this structure could be mapped 
onto a classic problem known to be in P (XOR-SAT). 

As we noted the Hamiltonians in the fc-local Pauli 
Commuting Hamiltonian problem can exhibit ground 
states that have topological entanglement [14| , An im- 
portant open question is whether this is true for other 
systems that can exhibit topological entanglement. The 
3-local Commuting Hamiltonian system doesn't possess 
topological order and is in NP [ll| . However 4- local Com- 
muting Hamiltonian systems do have topological order. 
Here the fc-local Pauli Commuting Hamiltonian system 
has topological order, but this problem was actually in P. 
We therefore conclude that the complexity of general fc- 
local Commuting Hamiltonian problem must depend on 
something more specific that whether or not a quantum 
system can exhibit topological order. 
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